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Abstract 

We discuss the effects of initial collision geometry and centrality bin definition on correlation 
and fluctuation observables in nucleus- nucleus collisions. We focus on the forward-backward cor- 
relation coefficient recently measured by the STAR Collaboration in Au+Au collisions at RHIC. 
Our study is carried out within two models: the Glauber Monte Carlo code with a 'toy' wounded 
nucleon model and the hadron-string dynamics (HSD) transport approach. We show that strong 
correlations can arise due to averaging over events in one centrality bin. We, furthermore, argue 
that a study of the dependence of correlations on the centrality bin definition as well as the bin size 
may distinguish between these 'trivial' correlations and correlations arising from 'new physics'. 

PACS numbers: 24.10.Lx, 24.60.Ky, 25.75.-q 
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I. INTRODUCTION 



Correlations of particles between different regions of rapidity have for a long time been 



considered to be a signature of new physics. A shortening in the corre 



ation length in rapidity 



Q. 



Conversely, the 



has been thought to signal a transition to a quark-gluon plasma 
appearance of long-range correlations has been associated with the onset of the percolation 
limit, also linked to the QCD phase transition 3|, |4j. Recently, the correlations across a 
arge distance in rapidity have also been suggested to arise from a color glass condensate 



,5, 



The observation of such correlations in A+A collisions at RHIC energies by the STAR 
Collaboration [7,8] has therefore elicited a lot of theoretical interest. 

The purpose of this work is to identify some baseline contributions to the experimentally 
observed correlations, contributions that do not depend on new physics. We will use models 
that incorporate event-by-event fluctuations in initial conditions: a 'toy' wounded nucleon 
model and the hadronic string dynamics (HSD) transport model to illustrate the effect 
of these contributions. We then argue that a study of the dependence of correlations on 
the centrality bin definition as well as the bin size may distinguish between these 'trivial' 
correlations and correlations arising from 'new physics'. 

The paper is organized as follows. In section II the main observables are introduced. In 
sections III and IV we study system size fluctuations and the resulting centrality dependence 
of correlations of two disconnected regions in momentum space within two different models: 
the Glauber Monte Carlo model (with no hadronic re-interactions or initial state dynamics) 
and the HSD transport model. Section V summarizes our study. 



II. DEFINITION OF OBSERVABLES 



The statistical properties of a particular sample of events can be characterized by a set of 
moments or cumulants of some observable. These properties depend upon a set of criterions 
which are used to select this sample. Applied to the context of heavy-ion collisions this 
translates to the construction of centrality bins of collision events from minimum-bias data. 
We will discuss the charged hadron multiplicities Na and Nb in two symmetric intervals Ar) 
of pseudo-rapidity. After construction of the centrality bins, one can calculate the moments 
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of a resulting distribution P^ 9ap (N A , N B ; Ar]): 

(N k A ■ N l B )*r> = £ N k A N l B P*»»(N A , N B] A v ) . (1) 

N A ,N B 

In Eq. ([T]) the subscript c denotes a particular centrality bin, while the superscript rj gap 
denotes the separation of two symmetric intervals Ar] in pseudo-rapidity space where particle 
multiplicities Na and N B are measured. The correlation coefficient 1 is defined by 

^ (AiV A ■ ANb)c 9 °" p 

^{{ANa) 2 )^ ((AN B )yr p 

and measures how strongly multiplicities Na and iVe - in a given centrality bin c for 
pseudo-rapidity separation r] gap - are correlated. In Eq. (J2D, AiV = iV — (N)^ sap and 
((AA^4) 2 )[? 9ap = ((AiV B ) 2 )2 sap for symmetric intervals. 

The recent preliminary data on forward-backward correlation coefficient (j2J) of charged 
particles by the STAR Collaboration [3, S| exhibit two striking features: a) an approximate 
independence on the width of the pseudo-rapidity gap r] gap , b) a strong increase of p with 
centrality. 



III. GLAUBER MONTE CARLO MODEL 

We use the PHOBOS Glauber Monte Carlo code [9j coupled to a 'toy' wounded nucleon 
model, referred to as GMC. The aim of this model is to emphasize two crucial aspects: 1) an 
averaging over different system sizes within one centrality bin introduces correlations; 2) the 
strength of these correlations depend on the criteria used for the centrality definition and 
on the size of the centrality bins. 

Employing the Glauber code we model the distribution of the number of participating 
nucleons, Np, in each nucleus-nucleus collision for given impact parameter b (cf. Fig. [TJ 
left). This is done for Au+Au with standard Wood-Saxon profile and the nucleon-nucleon 
cross section of <jnn = 42 mb. The 'event' construction proceeds then in a two-step process. 
Firstly, we randomly generate the total number of charged particles: 

N P 

N ch = E <h > (3) 
i=i 



1 We use a different notation from Refs. 0,0] denoting the correlation coefficient as p and reserve the letter 
b for the impact parameter. 



where the number of charged particles n l ch per participating nucleon are generated by in- 
dependently sampling a Poisson distribution with given mean value n c h = 10. Secondly, 
these charged particles are randomly distributed according to a Gaussian in pseudo-rapidity 
space: 

where ov, = 3 defines the width of the pseudo-rapidity distribution. Hence, in each single 
event there are no correlations between the momenta of any two particles. Note that nu- 
merical values of n c h and a v are fixed in a way to have a rough correspondence with the 
data on charged particle production at yfs = 200 GeV. 




FIG. 1: (Color online) Left: Histogram shows the distribution of events with fixed number of 
participating nucleons Np and fixed impact parameter b in Au+Au collisions at y/s = 200 GeV. 
Right: The scaled variance wp of the distribution of participating nucleons in 10% bins as defined 
via b, Np, and N^. 

In Fig. [U(Ze/£) we show the GMC event distribution in the (b, iVp)-plane. For each of these 
events we randomly generate the number of charge particles N c h and their redistribution 
according to Eqs. (j3J) and (j3J), respectively. The construction of centrality classes can now 
be done in several ways. Here we focus on the following criterions: via impact parameter 
b, via the number of participating (wounded) nucleons Np, and via the charged particle 
multiplicity N r jf in the midrapidity window \q\ < 1. 

In the case one chooses the number of participating nucleons Np for centrality definition, 
one takes vertical cuts in Fig. [T] (left), while choosing the impact parameter b, one takes 
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horizontal cuts. Hence, depending on the centrality definition, one may assign a particular 
event (characterized by Np and b) to two different centrality bins. 
In Fig. [T] (right) we show the resulting scaled variance u P , 



Up 



(Np)c 



(5) 



of the underlying distribution of the number of participating nucleons Np in each centrality 
bin. Using the centrality selection via impact parameter b, which is only the theoretically 
available trigger, one generally obtains a rather wide distribution of participating nucleons 
in each bin. The lines for centrality selections via and via Np are similar due to the 
event construction with Eqs. (I3ll4j) . An interesting feature of the GMC model is that ujp 
increases with centrality for the selection via Np. This conclusion of the GMC model seems 
to have a rather general origin. 

We now investigate the sensitivity of the forward-backward correlation signal as a function 
of the separation r] gap of two narrow intervals (A77 = 0.2) on the centrality definition. This 
is done for the 10% centrality defined via Np, via b, and via N r jf . The results are shown 
in Fig. [2 
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FIG. 2: (Color online) The forward-backward correlation coefficient p for 10% centrality classes 
defined via Np (left), via the impact parameter b (middle), and via the multiplicity in the central 
rapidity region N r jf (right). 



In the GMC we can identify the number of participating nucleons Np with the system 
size, and up as the measure for system size fluctuations. Having a large system as measured 
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by Np implies a large number of charged particles N c h- In GMC they are distributed 
independently in pseudo-rapidity space. Conversely, an event with small Np contains only 
few charged particles. By grouping the collision events into 10% centrality bins one finds 
rather large Ap-fluctuations in one specific bin. The averaging over different states in the 
centrality bin introduces correlations between any two regions of pseudo-rapidity. Small 
systems will have few particles 'on the left' and few particles 'on the right' with respect to 
midrapidity. Large systems will have many particles 'on the left' and many particles 'on the 
right'. But this just means a non-zero forward-backward correlations. From the definition 
(j2J) one finds a positive correlation coefficient p due to averaging over system sizes. 

Note that centrality selections via Np and via N c h give essentially the same results for p in 
the GMC (cf. left and right panels of Fig. [2]). Using the impact parameter b for the centrality 
definition generates centrality bins with almost constant p as seen in Fig. [2] (middle). This 
is due to a rather flat dependence of up on the centrality defined via b, as shown in Fig. [1] 
(right). In the GMC model the apparent ordering of p values with respect to centrality bins 
originates from the width of the underlying distribution in the number of wounded nucleons 
in each bin, i.e. from the values of up . 

The measured and apparently strong forward-backward correlations can be accounted for 
by a 'toy' model such as the GMC, provided it produces particles over the whole rapidity 
range and includes strong enough event-by-event fluctuations of Np. The next section will 
show that an introduction of dynamics and hadron re-interactions within the HSD does not 
alter these conclusions significantly. 



IV. HSD TRANSPORT MODEL SIMULATIONS 



A physically more reasonable scenario, which however also does not include any 'new 
physics' (such as color glass condensate, q uark-gluon plasma , et c . ) can be obtained in the 
Hadron- String-Dynamics (HSD) transport approach lOj, [ll|, ll2j . The HSD has been used 
for the description of pA, ttA and AA collisions from SIS to RHIC energies Q, fl. In 
this model, N, A, N*(1440), N*(1535), A, E and E* hyperons, 5, E* and Q as well as their 
antiparticles are included on the baryonic side, whereas the 0~ and l - octet states are in- 
corporated in the mesonic sector. Inelastic baryon-baryon (and meson-baryon) collisions 
with energies above y/s th ~ 2.6 GeV (and \fs th — 2.3 GeV) are described by the FRITIOF 
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string model 15j whereas low energy hadron-hadron collisions are modelled in line with 
experimental cross sections. As pre-hadronic degrees of freedom the HSD includes 'effec- 
tive' quarks (antiquarks) and diquarks (antidiquarks) which interact with cross sections in 
accordance with the constituent quark model. 




50 100 150 200 250 300 350 



FIG. 3: (Color online) The HSD and GMC distributions of events over Np. The vertical lines 
indicate 10% centrality bins. 



As before within GMC, the HSD events are generated according to a uniform distribution, 
N ev (b) ~ b. The resulting distribution of events in the (Np, 6)-plane is similar to the GMC 
result depicted in Fig. [1] (left). 

In Fig. [3] we show the distribution of events with fixed Np for both models. The vertical 
lines indicate 10% centrality bins as defined by the Np distribution. Note, that the peripheral 
part of the distribution determines also the centrality binning and the real bin widths. This 
is crucial for most central collisions where the number of events is small. Slight uncertainties 
in the peripheral "tail" of the distribution leads to large errors in the sizes of most central 
bins and hence to large changes in results for fluctuations and correlations. 

In contrast to the STAR data, we use in the HSD simulations the charged particle ref- 
erence multiplicity N r jJ in the same pseudo-rapidity range \q\ < 1 for all values of r\ gap . 
This procedure introduces a systematic bias, since the pseudo-rapidity regions for the mea- 
sured multiplicity in a small At] window (signal) and for the reference multiplicity partially 
overlap. This bias, however, is small and does not affect any of our conclusions. 

In Fig. H]we show the scaled variance of the underlying Np distribution for 10% (left) and 
2% (right) centrality bins defined via different centrality triggers within HSD. The results for 




FIG. 4: (Color online) The HSD results for the fluctuations Lu par t as a function of the mean value 
(Np) of the participating nucleons within bins as defined via b, Np, and N^. The left panel 
corresponds to a 10% and the right to a 2% bin width. 

10% bins can be compared with the scaled variance up in the GMC model in Fig. [1] {right). 
Fluctuations of the number of participants, as well as their average values, are similar in both 
HSD and GMC models when the centrality bins are defined via Np. These quantities are 
completely defined by the Np distribution, which is similar in both models (Fig. [3]). Binning 
via the impact parameter h in HSD, as well as in GMC, gives decreasing fluctuations in the 
participant number with increasing collision centrality. The results for 10% bins defined via 
the reference multiplicity are rather different in the GMC and HSD models. In GMC the 
charged multiplicity distribution is implemented according to Eqs. ( I3ll4"l) . Hence, the results 
obtained by binning via the reference multiplicity follow the line obtained by binning via Np. 
In contrast to the GMC, in the HSD simulations the average number of charged particles 
n c h per participating nucleon is not a constant, but increases with Np. Additionally, the 
shape of rapidity distribution is also different in different centrality bins. These two effects 
lead to different values of tup in the centrality bins defined via N r jf in the GMC and HSD 
models. 

One comment is appropriate here. It was argued in Ref. [if! that any centrality selection 
in nucleus-nucleus collisions is equivalent to the geometrical one via impact parameter b. 
This result was obtained in Ref. 16] by neglecting the fluctuations at a given value of b. 
Thus, different centrality selection criterions give indeed the same average values of physical 
observables. However, they may lead to rather different fluctuations of these observables in 
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the corresponding centrality bins, cf. equal values of (Np) and different values of uop for 
different centrality selections presented in Fig. HI 

When considering smaller centrality bins (2% in Fig. HJ right) the fluctuations in the 
participant number become smaller and more strongly dependent on the definition of the 
binning. 
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FIG. 5: (Color online) The HSD results for the forward-backward correlation coefficient p for 10% 
{top) and 2% {bottom) centrality classes defined via Np {left), via impact parameter b {center), and 
via the reference multiplicity N r l^ {right). The symbols in the top right panel present the STAR 
data in Au+Au collisions at ^/s = 200 GeV ?[ P 



The Fig. [5] summarizes the dependence of forward-backward correlation coefficient p as 
a function of r] gap on the bin size and centrality definition within the HSD model. The 
dependence of p on i] gap is almost flat, reflecting a boost-invariant distribution of particles 
created by string breaking in the HSD. The right top panel of Fig. [5] demonstrates also a 
comparison of the HSD results with the STAR data [7|, [§]. One observes that the HSD 
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results exceed systematically the STAR data. However, the main qualitative features of the 
STAR data - an approximate independence of the width of the pseudo-rapidity gap r\ gav and 
a strong increase of p with centrality - are fully reproduced by the HSD simulations. 

The correlation coefficient p largely follows the trend of the participant number fluc- 
tuations wp as a function of centrality. The actual results, however, strongly depend on 
the way of defining the centrality bins. For instance, choosing smaller centrality bins leads 
to weaker forward-backward correlations, a less pronounced centrality dependence, and a 
stronger dependence on the bin definition. The physical origin for this is demonstrated in 
Fig. El As the bin size becomes comparable to the width of the correlation band between 
Np and N^, the systematic deviations of different centrality selections become dominant: 
the same centrality bins defined by Np and by N r jf contain different events and may give 
rather different values of forward-backward correlations coefficient p. 




FIG. 6: (Color online) Histogram shows the distribution of HSD events with fixed number of par- 
ticipating nucleons Np and fixed reference charge particle multiplicity iV!? . The same centrality 
class (20-22% as an example) defined in various ways contains different events. 

It should be underlined that these properties are specific to the geometric nature of the 
correlations analyzed here. If the observed fluctuations are of dynamical origin (for exam- 
ple, arising from the quantum fluctuations of coherent fields created in the first fml c of the 

nn 

system's lifetime as in Refs. |5|,l6|), there are no evident reasons why they should strongly 
depend on centrality bin definitions and bin sizes. Thus, the experimental analysis for differ- 
ent bin sizes and centrality definitions - as performed here - may serve as a diagnostic tool 
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for an origin of the observed correlations. A strong specific dependence of the correlations 
on bin size and centrality definition would signify their geometrical origin. 



V. SUMMARY 

In conclusion, we have presented a study of the system size event-by-event fluctuations 
causing the rapidity forward-backward correlations in relativistic heavy-ion collisions. Our 
analysis has been based on two independent models - a 'toy' wounded nucleon model realized 
as a Glauber Monte Carlo event generator and the microscopic HSD transport approach. We 
have shown that strong forward-backward correlations arise due to an averaging over many 
different events that belong to one 10% centrality bin. In contrast to average multiplicities, 
the resulting fluctuations and correlations depend strongly on the specific centrality trigger. 
For example, the centrality selection via impact parameter b used in most theoretical cal- 
culations and via N r jJ used experimentally lead to rather different values of up and p and 
their dependence on centrality. 

In the HSD model the Np distribution is similar to that in the GMC It includes also 
the fluctuations in the number of strings and the fluctuations in the number of hadrons 
from individual string fragmentation. The HSD simulations reveal strong forward-backward 
correlations and reproduce the main qualitative features of the STAR data in A+A collisions 

fin 

at RHIC energies [5|, |6|] . 

The forward-backward correlations can be studied experimentally for smaller size central- 
ity bins defined by N r jf . When the size of the bins decreases, the contribution of 'geometri- 
cal' fluctuations discussed in our paper should lead to weaker forward-backward correlations 
and to a less pronounced centrality dependence. Let us stress that the 'geometrical' fluctu- 
ations discussed in our paper are in fact present in all dynamical models of nucleus-nucleus 
collisions. Thus, they should be carefully subtracted from the data before any discussion of 
new physical effects. We hope that a future experimental analysis in the direction exam- 
ined here will clarify wether the observed correlations by the STAR Collaboration at RHIC 
contain really additional contributions from 'new physics'. 
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